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We study the behaviour of five-dimensional fermions localized on branes, which we describe by do- 
main walls, when two parallel branes collide in a five-dimensional Minkowski background spacetime. 
We find that most fermions are localized on both branes as a whole even after collision. However, 
how much fermions are localized on which brane depends sensitively on the incident velocity and 
the coupling constants unless the fermions exist on both branes. 
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I. INTRODUCTION 



It has been known since the 70's that topological de- 
fects such as domain walls can trap fermions on their 
world volumes [lj. In the 80's this fact formed an in- 
tegral part of suggestions that one may regard our uni- 
verse as a domain wall @, [H, 0, [f| , or more generally a 
brane in a higher dimensional universe @, 0, B, Q ■ The 
idea is that the fermionic chiral matter making up the 
standard model is composed of such trapped zero modes 
EH El El El El EE El, G3- A similar mechanism 
is used in models, such as the Horawa-Witten model 
[3 El of heterotic M-Theory, in which two domain walls 
are present. Our world is localized on one brane and a 
shadow world is localized on the other brane. The exis- 
tence of models with more than one brane suggests that 
branes may collide, and it is natural to suppose that the 
Big Bang is associated with the collision [20l |2~H . This 
raises the fascinating questions of what happens to the 
localized fermions during such collisions? Put more pic- 
turesquely, what is the fate of the standard model dur- 
ing brane collision? In this paper we shall embark on 
what we believe is the first study of this question by 
solving numerically the Dirac equation for a fermions 
coupled via Yukawa interaction to a system of two col- 
liding domain walls, i.e. a kink-anti-kink collision in five- 
dimensional Minkowski spacetime. Each individual do- 
main wall may be described analytically by a static solu- 
tion and given such a solution one may easily find ana- 
lytically the fermion zero modes, which from the point 
of view of the 3+1 dimensional world volume behave 
like massless chiral fermions. The back reaction of the 
fermions on the domain wall is here, and throughout this 
paper, neglected. 

Kink-anti-kink collisions, have recently been studied 
numerically [H SI S3 . One solves the scalar field equa- 
tions with initial data corresponding to a superposition 
of the boosted profiles of a kink and an anti-kink. It 
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was found [22j,|23j that, depending on the initial relative 
velocity that such domain wall pairs can pass through 
one another, or bounce, or suffer a number of bounces 
in a fashion reminiscent of the cyclic universe scenario 
21 1 . One may extend the treatment to include gravity 
II [MS SI, SI] but in this paper we shall, for the sake 
of our preliminary study, work throughout with gravity 
switched off. One may now solve the Dirac equation in 
the time dependent background generated by the kink- 
anti-kink collision. We use as initial data for the Dirac 
equation the boosted profiles of the chiral zero modes 
associated with the individual domain walls. 

What we find was for us unexpected and quite remark- 
able. If the initial fermions exist on both branes, then 
without exception, for a whole range of initial conditions, 
the two initially distinct but localized fermion distribu- 
tions merge in the neighborhood of the collision, and then 
emerge after the collision again localized on one or the 
other kink. By contrast if one of the kinks is empty, which 
we refer to as a vacuum brane, then the amplitudes of the 
fermions on the kinks after the collision highly depend on 
the incident velocity and the coupling constants. 



II. FERMIONS ON MOVING BRANES 

A. Fermion with Yukawa coupling and its 
symmetry 

We start with a discussion of five-dimensional (5D) 
four-component fermions in a time-dependent domain 
wall in 5D Minkowski spacetime. As a domain wall, we 
adopt a 5D real scalar field $ with an appropriate po- 
tential V(Q). The 5D Dirac equation with a Yukawa 
coupling term gr$#^ is given by 

(T A d A +g*)9 = 0, (1 = 0,1,2,3,5), (2.1) 

where is a 5D four-component fermion. T A are the 
Dirac matrices in 5D Minkowski spacetime satisfying the 
anticommutation relations, 



{r A ,T B } = 2 V 



All 



(2.2) 
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where r\ AB — diag(— 1, 1, 1, 1, 1) is Minkowski metric 1 . 
We explicitly use the following Dirac-Pauli representation 



r° = 



-i 
i 

- 



p 5 = 







1 

1 

(k = 1,2,3) , 



(2.3) 



with <r k being the Pauli 2x2 matrices. 

Note that Eq. (|2.1|) implies current conservation law: 



d A n A = . 



(2.4) 



where n = vpr^vl/ is conserved number current. Here we 
define * = ^T 6 . This gives conserved number density 
n = n° = *r S * = * t *. The total number of fermions 
is defined by N = J d 5 Xn, which is conserved. 

Later we shall need the fact that the Dirac equation 
(|2.ip has the following time reversal and reflection sym- 
metries: 

(1) If ty(t,x,z) is a solution of the Dirac equation with 
scalar field x, z), T ^(—t, x, z) is a solution of the 
Dirac equation with the scalar field — $(— t, x, z), where 
X 5 = z is the coordinate of a fifth dimension. In par- 
ticular, when there is no interaction ($ = or g = 0), 
r ^(— t, x, z) is time reversal of ^(t, x, z) 

(2) If fy(t,x, z) is a solution of the Dirac equation with 
scalar field &(t,x, z), P 5 ^^, x, — z) is a solution of the 
Dirac equation with the scalar field — x, —z). In par- 
ticular, if x, z) is a solution for a kink [an anti-kink], 
T 5 ty(t, x, z) is a solution for an anti-kink [a kink]. It will 
turn out that the solution with a kink [an anti-kink] is 
related to positive [negative] chiral fermions, which are 
defined below (see next subsection). 

(3) Combining (1) and (2), we find that T h T a ^(-t, x, -z) 
is a solution of the Dirac equation with $(— t, x, —z) 

If we assume some symmetries for a domain wall, we 
find further properties for fermions as follows. 

(i) For the case of a static domain wall, (1) yields that 
T° i &(t,x',z) is a solution for an anti-kink [a kink] if 
ty(t, x, z) is a solution for a kink [an anti-kink]. 

(ii) If a domain wall is described by a kink (or an 
anti-kink), which has symmetry such that x, — z) = 
—<&{t,x,z), (2) yields that T 5 ^(t, x, z) is a solution for 
an anti-kink [a kink] if ^(t, x, z) is a solution for a kink 
[an anti-kink]. 

(iii) We may also have time symmetry such that 
£, x, z) = $(t,x,z) for collision of two walls. In fact 

we find from numerical analysis that this ansatz is ap- 
proximately correct [23T |. Assuming ^-reflection symme- 
try as well, we find from (3) that r S r 6 *±(-i, x, -z), 
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which is time reversal and z-reflection of "i>±(t,x, z), is 
also a solution for the same scalar field $(i, x, z). 

Before going to analyze concrete examples, we intro- 
duce two chiral fcrmion states 



i(i±r«)* 



(2.5) 



This definition implies 



I(l ± P*)* 



±=*±, \ (ITH *±=0. (2.6) 



Using the representation (|2.3p . we have 



where ip+ and tp_ are two-component spinors. 
The Dirac equation (|2.1[) is now reduced to 

(±d 6 +g9)^ ± +T%il> T = 0. 



B. Fermions on a kink (or an anti-kink) 



(2.7) 



(2. 



As for a domain wall, now we assume the potential 
form is given by V(<fr) = j ( < I >2 — r; 2 ) . Here we recall 
the dimension of some variables. Since we discuss five 
dimensional spacetime, we have the following dimension- 
ality: 

[4] = [ry] = L- 3 / 2 , [*] = L~ 2 , [g] = L l '\ [A] = L , (2.9) 

where L is a scale length. In what follows, we use units 
in which m I? (= ry 2 / 3 ) = 1. 

Then a domain wall solution is given by 

$ 



etanh ( — ) 



(2.10) 



where e = ± correspond to a kink and an anti-kink solu- 
tions and D = \j2j\ is the width of a domain wall. Note 
that &(z) is an odd function of z. 

As for a fermion, in the case of a static domain wall, 

separating variables as ip+ = ip + (x A ')/+(z) and ^_ = 

(4) 

ip -(a^ M ) / '—{z) and assuming massless chiral fermions on 

a branc, i.e. T^d^ip ±(x fJ ') = 0, we find the equations for 
f±(z) as 



(±%+5*(«))/± = 
With Eq. (|2.10[) . we find the solutions are 



f± « 



cosh ( -^J 



(2.11) 



(2.12) 



Note that the fermion wave function is an even function of 
z. Hence the positive-chiral (the negative-chiral) fermion 
is localized for a kink (an anti-kink ) but is not localized 
for an anti-kink (a kink). 
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To fix numbers of fermions on a wall, f± should be 
normalized up to an arbitrary phase factor </>±(o), which 
is set to be zero. Using a number density of fermions 
given by 

n = = 2 + tplip-) , (2.13) 

we normalize the total number of fermions localized on 
a static domain wall to be unity, i.e. N = 1. More 
precisely, for a kink (an anti-kink), we impose 



which gives 



f±(z) 



r( 9 D + ±) 



2^DT{gD) 



n±dz 


= 1, 


(2.14) 


1/2 




-gD 




cosh ( — ] 
\DJ. 


. (2.15) 



Using this solution, we can describe the wave function 
of fermion localized on a kink (or an anti-kink) as 



■^ + {x)f+{z) 



^ A \x,z)=\ t 4 7 (x)/ - (z) 
-t/j_(x)f^(z) 



(2.17) 



To quantize the fermion fields, we define annihilation 
operators of localized fermions on a kink and on an anti- 
kink by 

a K = (¥ K \V) and a A = (* (A) ,*) (2.18) 

Note that those two states are orthogonal, i.e. 
(3r(K)^(A)} = Q 



C. Fermion wave function on a moving domain wall 

To discuss fermions at collision of branes, we first dis- 
cuss fermions on a domain wall moving with a constant 
velocity. When a domain wall is moving, however, $ is 
time-dependent, and then the above prescription (sepa- 
ration of the fifth coordinate) to find wave functions is 
no longer valid. 

Since 3-space is flat, we expand the wave functions by 
Fourier series as 



1 



d A k e lkx ip±(t,z;k) . (2.19) 



(2tt) 3 / 2 

We find the Dirac equations become 

{±d 5 + ff $) j> ± ~ (id Q ± (fc • a)) V T = ■ (2.20) 

In what follows, we shall consider only low energy 
fermions, that is, we assume that k w 0, that is \k\ is 



enough small compared with the mass scale of 5D fermion 
(g$). The equations we have to solve are now 



id ip± = (t<9 5 + 5$) 1>i 



(2.21) 



Since up- and down-components of ip± are decoupled, 
we discuss only up-components here. Note that tak- 
ing into account k mixes the up- and down-components. 
With this ansatz, we can describe fermion by two single- 
component chiral wave functions as 



ip+{z,t) + 




r/>-{z,t). (2.22) 



For a localized fermion on a static kink (or an anti-kink) , 
the wave functions are ip±(z,t) — f±{z). 

Next we construct a localized fermion wave function 
on a moving domain wall with a constant velocity v. In 
this case, we can find the analytic solution by a Lorentz 
boost. We find for a kink with velocity v, 

4,f ] (z, t; v) = \[^~¥ K) h(z - vt)) , (2.23) 



4, {K) (z,t-v) 



and for an anti-kink with velocity v, 

£ A \z,t-v) = ^^(T^-vi)), 



ip ( +\z,t;v) 



2 

*7+ 1 



(2.24) 



1±1^) ( 7 (z - vt)) 



where ^{z) = f + (z) and ^{S) = f_(z) are static 
wave functions of chiral fermions localized on static kink 
and anti-kink, respectively, and 7 = — v 2 is the 

Lorentz factor. 

We can check that the total number of fermions is pre- 
served also in the boosted Lorentz frame. From Eqs. 
(|2.23[) and (|2.24|) . we find that n — 7 n. Integrating it in 
the 2-direction, we find 



dz n = I dz (7(2 — vt)) 

t— const J t— const 



dz n (z) = 1 . 



(2.25) 



If a domain wall is given by a kink [an anti-kink], 
we have only the positive-chiral fermions in a comov- 
ing frame [the negative-chiral fermions]. However, from 
Eqs (|2.23j) and (|2.24|) . we find that the negative-chiral 
modes [positive-chiral modes] also appear in this boosted 
Lorentz frame. For a kink, the ratio of number density 
of the negative-chiral modes to that of the positive-chiral 
ones is given by 7 2 w 2 / (7 + l) 2 . 

The above wave functions on a moving domain wall 
with constant velocity can be used for setting the initial 
data for colliding domain walls. 
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III. FERMIONS ON COLLIDING DOMAIN 
WALL 

A. Initial setup 

We construct our initial data as follows. Provide a kink 
solution at z — — zq and an anti-kink solution at z = Zq, 
which are separated by a large distance and approaching 
each other with the same speed v. We can set up as an 
initial profile for the scalar field 

<f>(z,t) = $ (K) (z,t;w) + <5> {A \z,t;-v) - 1, (3.1) 

where 

$ (K ' A) (z, t;v) = ± tanh( 7 (z - vt) /D) (3.2) 

are the Lorentz boosted kink and anti-kink solutions, re- 
spectively. Here we have chosen that the initial time is 
t = tin = —zq/v. The domain walls collide at t = 0. 

For fermions on moving walls, we first expand the wave 
function as 

* = (x, z; v)a K + ^ (x, z; -v)a A + 9™ (x, z)a B , 

(3.3) 

where (x, z; v) and (x, z; —v) are the wave func- 
tion of right-moving localized fermion on a kink and those 
of left-moving one on an anti-kink, respectively, which 
are explicitly by Eq. (j!T2"3"]) and Eq. ([272"1]) . We also 

denote the bulk fermions symbolically by W in ' (x, z). We 
do not give its explicit form because it does not play 
any important role in the present situation. We have 
assumed in Eq. that {^\x,z;v),9^ ) (x,z;-v) 

(B) 

and ty\ n (x, z)} form a complete orthogonal system. Note 

that {9^\x,z; v) and 9^\x,z;—v)} are orthogonal. 

Now we can set up an initial state for fermion by 
creation-annihilation operators. We shall call a domain 
wall associated with fermions a fermion wall, and a do- 
main wall in vacuum a vacuum wall. We shall discuss 
two cases: one is collision of two fermion walls, and the 
other is collision of fermion and vacuum walls. For initial 
state of fermions, we consider two states; 

|KA) = <44|0) (3.4) 
|K0) = 4|0) (3.5) 

where |0) is a fermion vacuum state. 

B. Outgoing states and expectation values 

We discuss behaviour of fermions at collision. After 
collision of two domain walls, each wall will recede to in- 
finity with almost the same velocity as the initial one 
v. Therefore we expect that positive chiral fermions 
stay on a left-moving kink and negative ones on a right- 
moving anti-kink. Those wave functions are given by 



V&out (x, z\ — v) and 9 on l(x, z;v). There may be bulk 
fermions which are left behind after collision, which wave 

( B) 

function is symbolically written by 9 ou ^(x,z). Since the 
initial wave functions (^^\x,z;v) and 9^\x, z; —v)) 
are described by the finial wave functions (*out ( x > z 'i ~ v ) > 
^i n l(x, z; v), and ^ut^: z) at t = t out = zq/v), we find 
the relations between them by solving the Dirac equation 
(12 . 2 If) . Those relations can be written as 

*£\x,z;v) ~ a K ^l(x,z;-v) + K ¥^l(x,z;v) 

+7k*£t (*,*). (3.6) 
9£\ X) z;-v) ~ a A 9itl(x,z;v)+p A ^l(x,z;-v) 
+ lA ^l(x,z). (3.7) 

In order to define final fermion states, we also describe 
the wave function as 

* - 9%l(x, z; -v)b K + 9£l(x, z; v)b A + *W(x, z)b B , 

(3.8) 

where &k, and 6b are annihilation operators of those 
fermion states. From Eqs. f3l3]) . 1(3^6]) . ([377) and (|3~8)) . 
we find 

b K = a K a K + f3 A a A , (3.9) 
6 A = "AdA + /3kOk (3.10) 

Using the Bogoliubov coefficients «k, /3k and a A ,(3 A , 
we obtain the expectation values of fermion number on 
a kink and an anti-kink after collision as 

(N K ) = (KA|6j,6 K |KA) = |« K | 2 + \(3 A \ 2 (3.11) 
(iV A > = (KA|6^6 A |KA) = |« A | 2 + |/3 K | 2 (3.12) 

for the case of |KA). If the initial state is |K0), we find 

(JV K ) = (K0|6^6 K |K0) = |a K | 2 (3.13) 
(N A ) ee (K0|6 A 6 A |K0) = |/3 K | 2 . (3.14) 

C. Time evolution of fermion wave functions 

In order to obtain the Bogoliubov coefficients, we have 
to solve the equations for domain wall $[23] and fermion 
9 numerically. For the time evolution of 9, we use the 
Crank-Nicholson method since it is generally shown to be 
useful for the parabolic type of partial differential equa- 
tion. 

In our simulation of two-wall collision, we have three 
unfixed parameters, i.e. a wall thickness (D) and an ini- 
tial wall velocity (v) and a coupling between fermions 
and a domain wall (g). From the solution (|2.15[) . we find 
the fermions are localized within the domain wall width 
D if g > 2/D. When g < 2/D, fermions leak out from 
the domain wall. Hence, in this paper, we analyze for the 
case of g > 2. We set D = 1, but leave v free. 



5 



Before showing our results for fermions, we summarize 
the behaviours of domain walls discussed in [23| . We 
find a bounce or a few bounces at the collision of domain 
walls, which depends in a complicated way on the initial 
velocityfThere is a fractal structure in the initial velocity 
space |22|). After the collision, two domain walls recede 
into infinity with almost same velocity ±u. It is similar 
to collision of solitons. 

To obtain the Bogoliubov coefficients, we solve the 
Dirac equation for the collision of fcrmion- vacuum walls, 
i.e. fermions arc initially localized on one wall, and the 
other wall is empty (^[^ (x, z; v) or ^\^'{x, z; —v) ). 

We shall give numerical results only for the case that 
positive chiral fermions are initially localized on a kink 
{9^\x,z;v)). Because of z-reflection symmetry dis- 
cussed in S III A[ we find the same Bogoliubov coefficients 
for the case that negative chiral fermions are initially 

localized on an anti-kink (^^\x,z;—v)), i.e. |q;k| 2 = 
|a A | 2 and |/3 K | 2 = |/?a| 2 . 

Setting g = 2 and v = 0.8, we show the result in Fig. 
[TJ The other chiral mode appears at collision and the 
wave function splits into two parts after collision. 
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(a) t = (initial) 
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(b) t — 32 (at collision) 
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(c) t = 56 (after collision) 

FIG. 1: Snap shots of the number density of the wave func- 
tion (n = \f r '\f r ) and those of two chiral states n± for collision 
of fermion- vacuum walls. We set D — 1, g = 2, and v = 0.8. 



the initial value and the final one (|7k| 2 = 1 — (I^kI 2 + 
|/3k| 2 )) corresponds to the amount of bulk fermions left 
behind. 

The Bogoliubov coefficients depend on the initial wall 
velocity. In Table U we summarize our results for differ- 
ent values of velocity. 



V 


9 = 2 


9 = 2.5 


|qk| 2 |/3k| 2 | 7 k| 2 


|a K | 2 \/3k\ 2 | 7 k| 2 


0.3 


0.94 0.056 0.004 


0.47 0.53 0.00 


0.4 


0.87 0.12 0.01 


0.57 0.40 0.03 


0.6 


0.69 0.30 0.01 


0.78 0.17 0.05 


0.8 


0.42 0.55 0.03 


0.88 0.02 0.10 



TABLE I: The Bogoliubov coefficients of fermion wave func- 
tions localized on each domain wall after collision ( |aK| 2 and 
|/3k| 2 ) with respect to the initial velocity v. We also show 
the amount of fermions escaped into bulk space (|7k| 2 = 



(|a K | 2 + |/3 K | 2 )). 



= 2.5 in Table H For 
2 and |/?k| 2 are almost 



We also show the case of g - 
the coupling constant g = 2, \cxk 
equal (0.44 and 0.55), but for g = 2.5, most fermions 
remain on the kink (|a K | 2 = 0.88 and |/3 K | 2 = 0.02). We 
find that the Bogoliubov coefficients depend sensitively 
on the coupling constant g as well as the velocity v. In 
Fig. [21 we shows the ^-dependence. 
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FIG. 2: The Bogoliubov coefficients (|a K | 2 , |/3k| 2 ) with v = 
0.4 in terms of a coupling constant g. The circle and the 
cross denote |ok| 2 and |/3k| 2 respectively. Two sine curves 
(|a K | 2 , |/3k| 2 w [l±sin(4.2g-1.2))]/2 show the formula (1338]) 
with the best-fit parameters. 



From the asymptotic behaviour of the wave function as 
t —> oo, we obtain the Bogoliubov coefficients numerically 
such that |ok| 2 = 0.42 and \Pk\ 2 = 0.55. Since a few 
amount of fermions escapes into bulk space at collision, 
|o:pc| 2 + |/3k| 2 is not conserved, and the difference between 



Since the wave function is changed at collision, when 
the background scalar field evolves in a complicated way, 
one might think that the behaviour of wave function 
would be difficult to describe analytically. However, we 
may understand the qualitative behaviour in terms of the 
following naive discussion. 
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Before collision, the wave function is approximated 
well by ^\x, z; v). In order to evaluate the wave func- 
tion of fermion after collision, we have to integrate the 
Dirac equation (|2.21[) . During the collision, the spatial 
distributions of fermion wave functions are well-described 
by some symmetric function of the z-coordinate (see Fig. 
Q] (b)). So we may approximate them as 



V±=^±(*)e #±( *Vo(z), 



(3.15) 



where ipo( z ) is a normalized even real function. A± and 
(j)± are regarded as the amplitudes of positive- (negative- 
) chiral modes and those phases, respectively. The scalar 
field $ evolves as $ : 1 — > $ c (« —1.5) — > 1 at the col- 
lision point [z = 0). If we approximate the scalar field 
as 4> = $ c at collision for collision time Ai(~ D/c), in- 
tegration of Eq. (|2.21[) with respect to z gives the change 
of amplitudes and phases of wave functions as 



1 



1-A% 



d (f>± = -g$ c 



-.d A± = ±g<f> c sm(A<j>) , (3.16) 



\-A\ 



A± 



■cos(A0), (3.17) 



where A(j> = <j>_ — <f)+ . We have also assumed that total 
amplitude of wave functions is normalized = 1). 

This means that we ignore bulk fermions, which may be 
justified because |7k| 2 <C 1. If A<f) = and A± = 1 
initially, then we find A + (At) = 1 (or A_(At) — 1) from 
Eq. (|3.16p . which guarantees A<p = anytime from Eq. 

(SHU). 

We find that {A+,A(j>) = (1,0) (or (A_, A</>) = (1,0 )) 
is a fixed point of the system (Eqs. (|3. 16[) and (|3.17l) '). 
However, it turns out that those are unstable. On the 
other hand, we find that Acj) = ir/2 (or — 7r/2) is an 
attractor (stable fixed points) of the present system. 
The time scale to approach these attractors is given by 
(ff|* c |) _1 if AA 2 = A 2 _ - A\ = 0(1). 

Once we assume A</> = ±7r/2, then we find that the 
phases cj>± do not change. Then we can integrate Eq. 
([OB]) , finding 



1 



A ± (At) = - [1 ± sin {2eg<$> c At + C )] 



(3.18) 



where e — ±1 and Co is an integration constant. 

This formula may provide a rough evaluation of 
|o<k| 2 ) I /3k 1 2 • Comparing the numerical data and the for- 
mula p,18p with <f> c w —1.5, we find the fitting curves 
in Fig. H (e = -1, At w 1.4 and C = -1.2). The 
above naive analysis explains our results very well. We 
then conclude that A(j> = ±ir/2 is generic except for a 
highly symmetric and fine-tuned initial setting (A + = 1 
or A_ = 1 and Acj) — 0), and the formula (|3 . 18[) with 
A<p = ±7r/2 is eventually found after collision. The small 
difference may be understood by the details of the com- 
plicated dynamics of colliding walls. 



D. Fermion numbers on domain walls after collision 

We can evaluate the expectation values of fermion 
numbers after collision as follows. For the initial state 
of fermions, we consider two cases: case (a) collision of 
two fermion walls |KA) and case (b) collision of fermion 
and vacuum walls |K0). 

In the case (a), we find 

(N K ) = \a K \ 2 + \/3 A \ 2 = \a K \ 2 + \/3 K \ 2 

= 1-| 7K | 2 «1 (3.19) 

(N A ) = \a A \ 2 + |/? K | 2 = |«a| 2 + |/3a| 2 

= 1-| 7 a| 2 «1. (3.20) 

We find that most fermions on domain walls remain on 
both walls even after the collision. A small amount of 
fermions escapes into the bulk spacetime at collision. 
In the case (b), however, we obtain 



<^k) = |«k| 2 , (N A ) = |/? K | : 



(3.21) 



Since the Bogoliubov coefficients depend sensitively on 
both the velocity v and the coupling constant g, the 
amount of fermions on each wall is determined by the 
fundamental model as well as the details of the collision 
of the domain walls. 



IV. CONCLUDING REMARKS 

We have studied the behaviour of five-dimensional 
fermions localized on domain walls, when two parallel 
walls collide in five-dimensional Minkowski background 
spacetime. We have analyzed the dynamical behavior of 
fermions during collision of fermion-fermion branes (case 
(a)) and that of fermion- vacuum ones (case (b)). 

In order to evaluate expectation values of fermion num- 
ber on a kink and an antikink after collision, we solve the 
Dirac equation for the wave function in the case (b) and 
find the Bogoliubov coefficients, in which /3k denotes the 
amount of fermions transfering from a kink to an an- 
tikink (a vacuum wall). As a result, in the case (b) some 
fermions jump up to the vacuum brane at collision. The 
amount of fermions localized on which brane depends 
sensively on the incident velocity and the coupling con- 
stants g/X where g and A are the Yukawa coupling con- 
stant and that of the double-well potential, respectively. 
It can be intuitively understood that the amount of lo- 
calized fermions is roughly determined by the duration of 
collision for which they transfer to another wall or stay 
on the initial wall, and the localization condition depend- 
ing on the Yukawa coupling constant between fermions 
and domain walls. 

On the other hands, in the case (a), we find that most 
fermions seem to stay on both branes even after collision. 
This is because of the relationship \/3k\ = \Pa\, which is 
guaranteed by a left-right symmetry in the present sys- 
tem. This result means physically that the same state 
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(k ~ 0) of fermions are exchanged for each other by the 
same amount. Therefore, the final amounts of fermions 
does not depend on parameters. 

We conclude with some comments about the subject 
not mentioned above: 

(1) For the case of g < 2/D, the localization of fermions 
on a domain wall is not sufficient. The tail of fermion 
distribution extends outside the wall. As a result, we 
find that a considerable amount of fermions escapes 
into a bulk space at collision. For example, we find 
|a K | 2 + |/?k| 2 = 0.64 for g = 1 and v = 0.8. The for- 
mula (|3.18p is also no longer valid in this case (see Fig. 
[2]). This is because localization is not sufficient. 

(2) The collision of domain walls is rather complicated. 
We find a few bounces at collision depending on the inci- 
dent velocity. The number of bounces is determined in a 
complicated way (a fractal structure in the initial phase 
space [111, [23j]). Thus if we change the incident velocity 
very little, the number of bounces changes. This causes 
a drastic change of final distribution of fermions on each 
wall for the case (b). 

(3) Since we have discussed only the case of zero- 
momentum fermion on brancs (fc = 0), we have only a 
single state on each brane, which constrains the fermion 
number to be less than unity. If we take into account 
degree of freedom of low energy fermions, we can put dif- 
ferent states of fermions on each brane. As the result, 
the final state of fermions after collision is different from 
the initial state, and it depends sensitively on the cou- 
pling constant as well as the initial wall velocity just as 
the case of collision of fermion- vacuum walls. 

(4) In the case of collision of two vacuum branes, nothing 
happens in the present approximation. The pair produc- 



tion of fermion and antifermion, for which we have to 
take into account the momentum k, may occur at colli- 
sion. This pair production process may aslo be important 
in the cases of collision of two fermion branes and that 
of fermion- vacuum branes. The work is in process. 

(5) Inclusion of self-gravity is important. It changes the 
fate of domain wall collision [25J] . It would be interesting 
to see what happens to the fermion distribution when 
we have a singularity. Although we are now analyzing 
it based on a supergravity model [3p| . it may be more 
important to study the model based on superstring or 
M-theory. 

We will publish the results elsewhere. 
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